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Abstract 

Based on the unified analyses of KASCADE, AKENO, EAS-TOP and ANI EAS size spec- 
tra, the approximations of energy spectra of different primary nuclei have been found. The 
calculations were carried out using the SIBYLL and QGSJET interaction models in 0.1-100 
PeV primary energy range. The results point to existence of both rigidity-dependent steepen- 
ing energy spectra at R ~ 200 — 400 TV and an additional proton (neutron) component with 
differential energy spectrum (6.1 ± 0.7) ■ 10 -11 (E/Et)- 1 - 5 (m 2 -s-sr-TeV) _1 before the knee 
E k = 2030 ± 130 TeV and with power index 72 = -3.1 ± 0.05 after the knee. 
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1 Introduction 



High statistical accuracy in modern EAS experiments in the knee region encouraged the investi- 
gation of the fine structure of EAS size spectra. Although the origin of the knee is still a matter 
of debate, recently the series of publications Q appeared, where the sharpness and the spectral 
structure in the knee region were interpreted by the contribution of heavy nuclei from a single 
local supernova. Along with this, the absolute differential EAS size spectra measured at different 
atmosphere depths and different zenith angles are not explained yet from the point of view of 
a single A — Aai T interaction model and a single model of primary energy spectra and elemen- 
tal composition. Such an attempt has been made in Q based on the QGS interaction model 
and rigidity-dependent steepening primary energy spectra |j| for a description of vertical MSU, 
AKENO and Tien-Shan EAS size spectra. In the present work we worked out a formalism of the 
inverse problem solution - reconstruction of the primary energy spectrum and elemental composi- 
tion based on the known EAS size spectra of KASCADE {§, AKENO §], EAS-TOP j| and ANI 
[0 measured at different zenith angles. The calculations were done in the frames of QGSJET Q 
and SIBYLL j^] interaction models. As a primary spectrum we have tested the modified rigidity- 
dependent steepening primary energy spectra and the hypothesis of the additional component in 
the knee region jnj, O, [l^, [l3| . In this case the type of nucleus of the additional component was 
considering as unknown and determining by the best fit of the fine structure of EAS size spectra 
in the knee region. 



2 EAS inverse problem 

In general, the energy spectra (d^A/dEo) of primary nuclei (A) and detectable EAS size spectra 
(dl/dN*) are related by the integral equation - 

where Eq,A,9 are energy, nucleon number (1-59) and average zenith angle of primary nuclei, E e 
is an energy threshold of detected EAS electrons, N*(E > E e ) is the estimation value of EAS 
size obtained by the electron lateral distribution function. Here, by the EAS size (N e (E > 0)) we 
mean the total number of EAS electrons at given observation level (t). The kern (We) of integral 
equation (1) is determined as 

1 f 02 r dn(E ,A,e,t) d*(N e ) . ...... 

We = 7T / <w 6 7r\rZ 9d8dN e 

A J 9l J Q dN e dN* 

where dVl/dN e is an EAS size spectrum at the observation level (t) for given Eq,A,9 parameters 
of a primary nucleus and depends on A — Aait interaction model; A# = cos Q\ — cos 9i ; 



P 9 = P(Ne,E ,A,8) = ^-y J J D(N e ,E a ,A,9,x,y)dxdy 



is a probability to detect an EAS by scintillation array at EAS core coordinates \x\ < A/2, 
\y\ < Y/2 and to obtain the estimations of EAS parameters (A*, s - shower age, x*,y* - shower 
core location) with given accuracies; /dN* is a distribution of N*(N e , s, x, y) for given EAS size 

(Ne)- 

One may achieve significant simplification of equation (1) providing the following conditions during 
experiments: 

a) selection of EAS cores in a range where Pg = 1, 

b) the log-Gaussian form of the measuring error (d^ /dN*) with an average value ln(A e • 5) and a 
RMSD (Tat, where 5 involves all transfer factors (an energy threshold of detected EAS electrons, 7 
and ji contributions) and slightly depends on E and A, 

c) transformation (standardization) of the measured EAS size spectra to the EAS size spectra at 
observation level 

dl(0,9,t) „ dI(E e ,9,t) 
dN e ~ V dN* ' 

where r\ = f^ 7 " -1 ' cxp{(7 e — l) 2 a^/2} and 7 e is the EAS size power index, 

d) consideration of either all-particle primary energy spectrum d^sz/dE with effective nucleus 
A e ff(E ) or energy spectra of primary nuclei (d^SA^/dEo) gathered in a limited number of groups 
(£ = 1, . . . £ max ) as unknown functions. 

(a-d) conditions make EAS data of different experiments more comparable and equation (1) con- 
verts to the form 



di(o,e,t) _ f°° d^dn(E a ,A eff (E ),e,t) 



poo 

V \ 

JE mi 



or 



dN e 'J~ 8E 8N e 



dE (2) 



dN e r, ^J Emin dE 8N e 
However, even in this form the determination of primary energy spectra by measured EAS size 
spectra and solution of integral equations (2,3) in general is unsolvable problem. At the same time, 
using the a priori information about energy spectra of primary nuclei (d^A^/dEo) and the EAS 
size spectra dl/dN* = fi,j(N* f , 6j, t) measured in i = 1, . . . m size intervals and j = 1, . . . n zenith 
angular intervals, one may transform the inverse problem into ^-minimization problem 

-Mx->--{EE%3r-} <•« 

i j J 

with unknown (free) spectral parameters. Here Fij = F(N* i ,9j,t) are the expected EAS size 
spectra determined at the right hands of equations (1-3) and 07, of are the uncertainties (RMSD) 
of measured (fij) and expected (F iy j) shower size spectra. 
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One may also unify the data of different experiments applying minimization Xjj with re-normalized 
EAS size spectra 

where index k = 1, . . . I determines the observation levels (t) of experiments. Expression (5) offers 
an advantage for experiments where the values of methodical shift (5) and measuring error (ctat) 
are unknown or known with insufficient accuracy. 

The energy spectra of primary nuclei are preferable to determine (a priori) in the following gener- 
alized form 

Unknown (free) spectral parameters in approximation (6) are (3, Eknee(A) (so called "knee" of 
energy spectrum of A nucleus), 71 and 72 (spectral asymptotic slopes before and after knee), e 
(sharpness parameter of knee, f < e < 10). The values of $^4 and 71(A) parameters are known 
from approximations of balloon and satellite data |l4| at A = 1, 4, . . . 59 and Eq — 1 — 10 3 TeV. 
Parameter (3 ~ 1 determines the normalization of spectra (6) in 10 2 — 10 5 TeV energy range. 
Thus, minimizing \ 2 - functions (4,5) on the basis of measured values of dl(9i^k) / 'dN e j^ and corre- 
sponding expected EAS size spectra (2,3) at given m zenith angular intervals, n EAS size intervals 
and / experiments one may evaluate the parameters of the primary spectrum (6). Evidently, the 
accuracies of solutions for spectral parameters strongly depend on the number of measured intervals 
(m ■ n- 1), statistical errors and correctness of <9fi(i?o, A, 9, t)/dN e determination in the framework 
of a given interaction model. 



3 Results 

Here, the parametric solutions of the EAS inverse problem are obtained on the basis of KASCADE 
[§ (t =1020 g/cm 2 ), AKENO § (910 g/cm 2 ), EAS-TOP @ (810 g/cm 2 ) and ANI § (700 g/cm 2 ) 
published EAS size spectra. These experiments were carried out at different observation levels 
and were chosen for two reasons: satisfaction of (a-c) conditions from the section (2) and a high 
statistical accuracy of presented data (especially KASCADE experiment). Unfortunately, during 
the standardization of EAS size spectra (condition (c)) the value of r\ parameter is not always 
known with proper accuracy, which is the main reason of discrepancy in the results of different 
experiments. In our calculations we included r\ in the list of unknown spectral parameters and 
determined by the minimization of functional (4). The problem does not exist if the minimization 
of re-normalized EAS size spectra is unified, since the linear parameters (77^ • f3) are canceled out 
from the functional (5). 

The differential EAS size spectra dn(E , A, 9, t)/dN e for given E Q = 0.032, 0.1, 100 PeV, 
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A = 1,4,12,16,28,56, t = 0.5, 0.6, . . . , 1 Kg/cm 2 , cos<9 = 0.8,0.9,1 were calculated using COR- 
SIKA562(NKG) EAS simulation code |f§ at QGSJET [| and SIBYLL § interaction models. 
Intermediate values are calculated using 4-dimcnsional log- linear interpolations. The estimations 
of errors of the expected EAS size spectra dfl/dN e at fixed E , A, 9, t parameters did not exceed 
3 - 5%. 

The basic results of minimizations (4,5) at a given number (y) of unknown spectral parameters and 
the values of y 2 jq (or Xu/l u f° r un i ne d data), are presented in Tables 1-4, where q = mn — v — 1 
and q u — ^2 k (fnn)k — v — 1 are corresponding degrees of freedom. The upper (lower) rows of 
each experiment in Tables 1 and 4 correspond to parameters obtained by the QGSJET (SIBYLL) 
interaction model. 

3.1 Test of rigidity-dependent energy spectra 

Table 1 contains the approximation values of spectral parameters at rigidity-dependent approach 
1 

E knee {A) =R-Z (7) 

where Z is a charge of A nucleus and R is a parameter of magnetic rigidity (or a critical (cutoff) 
energy in more modern model |ll| 0). The results were obtained by minimizations (4,5) applying 
(6,7). The magnitudes of 71 (^4) for all nuclei were taken from Jl4| . The number of unknown (free) 
parameters is equal to v = 4 and corresponding solutions at two interaction models are presented 
in Table 1. 

The stability of solutions (or the steep of \ 2 global minimum) of minimizations (4,5) is seen from 
obtained errors. So, it is seen that the x 2 minimum for AKENO data does not depend on the e 
sharpness parameter at SIBYLL and partly at QGSJET interaction models. 

The obtained slopes of primary energy spectra after knee agree with the same calculations 
performed by QGS model and exceed well known expected values (3 — 3.1) in the ~ 10 17 eV energy 
range Such a steep of primary spectra after the knee is a result of (72, R) correlations in (4-7). 
In case of R ~ 2000 TeV the Ek(Fe) ~ 5.2 • 10 4 TeV and the primary energy spectra in the large 
interval [Ef.{H) — Ek(Fe)) at fixed 71 (Ag) can be conformed with corresponding EAS size spectra 
provided abnormally steep slope (72 — 3.4 — 3.5) after the knee. 

The results of expected EAS size spectra in comparison with corresponding experimental data are 
shown in Fig. 1 (the thin solid lines by QGSJET model, the thin dashed lines by SIBYLL model). 
Despite the satisfactory agreement (x 2 ~ 1) of EAS size spectra with predictions of rigidity- 
dependent steepening spectra (6,7) and QGSJET interaction model the form (fine structure) of 
the measured EAS size spectra in the knee region differs from the form of corresponding expected 
spectra. It is worth mentioning that the difference is formally small and does not exceed several 
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Table 1: Rigidity (R), slope (72), "sharpness" (e), shift (77 • /?) and corresponding x 2 /l values 
obtained by approximations of EAS size spectra at QGSJET (upper rows) and SIBYLL (lower 
rows) interaction models and rigidity-dependent assumption (7). 



Model 


E k (H)[TeV] 


E k (He,Li) 


E k (Be-Na) 


E k (Mg-Ci) 


E k (Ar-M) 
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Table 2: Spectral parameters E knee (A) [TeV] for different groups of primary nuclei and different 
interaction models. The results obtained by unified analyses of EAS data at 72 = 3.1 and e = 4. 

percents in precise KASCADE data. 

As a next step, we attempted to test the rigidity-dependent approximation (7) directly. The knee- 
energies Ek(A^) were chosen as unknown (free) parameters in the approximations of EAS data 
(4-7), where £ = 1, . . . £ max . However, the stable solutions for free parameters were obtained only 
at unified minimization (5), fixed values 71 (A^), 72,5 and number of nuclear groups £ max < 5. 
In Table 2 the values of spectral parameters Eknee(A^) obtained by minimization of Xu (5) for 
5 groups of primary nuclei (H), (He, Li), (Be — Na), (Mg — CI), (Ar — Ni) at given values of 
72 = 3.1 , e = 4.0 and v = 5 are presented. It is seen, that approach (7) is performed only for 
nuclei with A > 1 and R ~ 400 TV or (7) is valid for all nuclei at R ~ 400 TV but there is 
an additional proton flux with E^ ee ~ 3000 TeV which shifts the knee value of the total proton 
energy spectrum. 

The following testing of the rigidity-dependent approach (7) is based on the investigation of the 
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Tabic 3: Parameters of all-particle energy spectrum obtained by unified EAS data (the first and 
second rows). The last row corresponds to results of extrapolation of fits |l4[ ] taking into account 
the assumptions (6,7) at R = 600 TV. 



all-particle energy spectrum. Toward this end the fits of primary energy spectra d^SA/dEg for 
different nuclei (A = 1, . . . 59) known from [[uj were extrapolated up to Ea = 10 5 TeV energies 
taking into account (6,7) at R ~ 600 TeV |u|. The obtained expected all-particle energy spectrum 
dl-s/dEa — ^2 A d^SA/dEo was approximated by expression similar to (6) at five free parameters 
($£, 71, 72, Ek, e). The average values of primary nuclei (exp(lnA)) obtained from extrapolations 
of spectra [[li] were approximated by step function 

A eff (E ) = a + b-\n(^) (8) 

where 6 = b± at Eq < Ek and b = 62 at Eq > Ek ■ The results of these approximations are presented 
in a last row of Table 3 . The value of sharpness parameter was equal to e = 1 ± 0.1. 
The first and second rows of Table 3 content the parameters of the all-particle energy spectra 
(9Ss/9i?o), which were obtained by minimization Xu (expressions 2,5-7) of unified EAS size data 
at v — 6 and e = 1. The approximation (8) has been used for A e ff(Eo) which is at the right 
hand of expression (2). ft is necessary to note that the solutions for 0,61,62 parameters can be 
obtained only by re-normalized EAS size spectra (5) because of the strong correlation between a 
linear parameter 77 and an effective nucleus A e ff. 

It is seen from Table 3 that the model of rigidity-dependent energy spectra predicts the increase 
(bi > 0, third row of Table 3) of the mean nucleus with energy, whereas the presented analysis 
of EAS data points out a decrease of A e ff with energy (bi < 0, first two rows of Table 3) in the 
energy range of Eq < Ek- It is obvious, that despite A could not be exactly equal to A e ff, at least 
their dependence on energy must be the same. 



The results of recent precise experiments DICE 17 and C ASA-BLANK A |18| also point out to 
the decrease of In A with energy at Eq < Ek ■ This dependence of A e f f and In A on energy might 
be explained by the contribution of an additional light component in a primary nuclei flux. 



From the above analyses follows that rigidity-dependent steepening energy spectra in combina- 
tion with QGSJET or SIBYLL interaction model can not explain the obtained results of the fine 
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structure of EAS size spectra Q in the knee region (Table 1 and Fig. 1), the large values of 

knee Ek{H) for Hydrogen component (Table 2) and dependence of the effective nucleus A e ff(Ea) 
on primary energy before the knee (Table 3). In this connection we have carried on the search of 
more adequate model of primary energy spectra and elemental composition. 

3.2 Test of additional component 

|l3| the primary energy spectra in approximation (6) have been 
added by a new (polar cap component d^Add/dE® with power energy spectrum 




where = 7 } p) at E < E { k p > and 7W = 72 at E > Eg*. 

New spectral parameters <&( p \ j^ p \ E^ and nucleon number A&> of the additional component 
are considered as unknown and determined together with parameters of rigidity-dependent energy 
spectra (6,7) by minimization of x 2 and xfr (4,5) at v = 7. The results of expected EAS size spectra 
for each experiment (KASCADE, AKENO, EAS-TOP and ANI) taking into account contribution 
of additional component (9) are shown in Fig. 1 (the thick solid lines by QGSJET model and 
the thick dashed lines by SIBYLL model). It is seen that the additional component with high 
accuracy (2-5% for KASCADE and ANI data) describes the fine structure of EAS size spectra in 
the knee region. The values of slopes of additional component before the knee turned out to be 
7^ ~ 1.5±Q 2- The nucleon number (A^) of this component with high reliability did not exceed 
of A^ = 1 for most of experiments especially at QGSJET interaction model (except from AKENO 
(yJW ~ 56)). The unified analyses of all experiments at QGSJET and SIBYLL interaction models 
also gave a proton or neutron [Ay^ — 1) composition of the additional component. The values of 
other spectral parameters at 7^ = 1.5 and A^ = 1 are included in Table 4. 
The obtained result disagrees with [|l] where the alike component consists of several heavy nuclei. 
However, our result is based on the high accuracy of the coincidence (~ 2 — 3% for KASCADE 
data) of expected and measured EAS size spectra by both x 2 criterion and the overlapping of fine 
structure of spectra. 

The comparison of parameters 72 from Tables 1,4 shows that spectral slopes after the knee from 
Table 4 roughly overlap with the expected slope (3 — 3.1) well known from N e ^> 10 7 EAS data. 
This can be explained the fact that the spectral brake of summary proton component ~ 2 • 10 3 

TeV) is closer to the iron component {E k (Fe) ~ 10 4 TeV) at R ~ 400 TV. 

The coexistence of rigidity-dependent primary energy spectra and additional proton flux with 
spectral parameters E^ ~ 2000 TeV, 7^ ~ 1.5, 72 =72 —3.1 explains also the shifted value of 
the knee Ek(H) in Table 2 and the decrease of A e ff(Eo) at Eq < Ek (see section 3.1). 
Although the additional component does not contribute to shower sizes below or above the knee 



Based on predictions [ROl Ml, 
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Table 4: Spectral parameters taking into account the contribution of additional proton component. 



significantly , it is essential for the reproducing of the sharp knee of EAS size spectra. This idea is 
taken from ]l^] and is confirmed here. 

The final all-particle energy spectrum (dl/dEo) obtained by unified EAS size data at QGSJET 
interaction model 

dl _ t \ - d^ A d^Add \ 
dE Q dE dE Q J 

and corresponding energy spectra (d^$A/dEo) of 6 nuclear groups with additional component 
(d^$Add/dEo) at normalization (3 — 1 |Q are presented in Fig. 2. The solid (dashed) line is the 
all-particle energy spectrum obtained by unified (only KASCADE) EAS size spectra. Dotted lines 
are the energy spectra of different primary components obtained by unified EAS data. Symbols in 
Fig. 2 are the data from DICE 0, CASA-BLANKA §J] arrays and reviews H). 



The numerical values of all-particle energy spectrum (solid line in Fig. 2), corresponding error, 
spectrum of additional proton component and average nucleus versus primary energy are presented 
in Table 5. 



Conclusion 

High statistical accuracy of experiments KASCADE, EAS-TOP, AKENO and ANI allowed to 
obtain approximations of primary energy spectra and elemental composition with accuracy ~ 15% 
in the knee region. Along with this, KASCADE and ANI EAS size spectra at 6 < 37° are described 
with the accuracy of ~ 2 — 5% in a whole measurement interval. 
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0.163E-10 


0.33E-11 


0.252E-11 


2.05 


0.711E+04 


0.843E-11 


0.18E-11 


0.126E-11 


2.11 


0.888E+04 


0.435E-11 


0.10E-11 


0.634E-12 


2.16 


0.111E+05 


0.222E-11 


0.57E-12 


0.318E-12 


2.19 


0.217E+05 


0.280E-12 


0.77E-13 


0.402E-13 


2.19 


0.827E+05 


0.448E-14 


0.14E-14 


0.642E-15 


2.19 



Table 5: Expected all-particle flux, error, flux of additional component and average nucleus 
different primary energies. 
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Obtained results show the evidence of QGSJET interaction model at least in 10 5 — 10 7 TeV energy 
range, rigidity-dependent steepening primary energy spectra at R ~ 200 — 400 TV and existence 
of the additional proton (or neutron) component with spectral power index 7 < p) ~ 1.5±g;| before 
the knee E^ ~ 2030 TeV. The contribution of the additional proton (neutron) component in 
all-particle energy spectrum turned out to be 20 ± 5% at primary energy E = E^\ 
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EAS size (N.) 

Figure 1: KASCADE, AKENO, EAS-TOP and ANI EAS size spectra (symbols). The thin (thick) 
lines correspond to predictions via rigidity-dependent steepening primary spectra (with the addi- 
tional proton component). The solid and dashed lines correspond to the QGSJET and SIBYLL 
interaction models respectively. 



13 




c / 

, A.c. 



7 3 4 

10 10 10 

E (TeV/ particle) 



Figure 2: Primary energy spectra and elemental composition. The solid (dashed) line is the all- 
particle energy spectrum obtained by unified (only KASCADE) EAS data. The dotted lines are 
the energy spectra of different nuclear groups. The A.c. dotted line is the energy spectrum of 
additional proton component. The symbols are the data from |l7], [l8| and reviews [^4[ fl6| . 
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